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ABSTRACT 

We present calculations of the heliospheric solar wind charge exchange (SWCX) emission spectra 
and the resulting contributions of this diffuse background in the ROSAT j keV bands. We compare 
our results with the soft X-ray diffuse background (S XRB) emission detec ted in front of 378 identified 
shadowing regions during the ROSAT All-Sky Survey ( Snowden et al.l200ol) . This foreground component 
is principally attributed to the hot gas of the so-called Local Bubble (LB), an irregularly shaped cavity of 
~ 50-150 pc around the Sun, which is supposed to contain ~10 6 K plasma. Our results suggest that the 
SWCX emission from the heliosphere is bright enough to account for most of the foreground emission 
towards the majority of low galactic latitude directions, where the LB is the least extended. On the 
other hand, in a large part of directions with galactic latitude above 30 degrees the heliospheric SWCX 
intensity is significantly smaller than the measured one. However, the SWCX R2/R1 band ratio differs 
slightly from the data in the galactic center direction, and more significantly in the galactic anti-centre 
direction where the observed ratio is the smallest. Assuming that both SWCX and hot gas emission are 
present and their relative contributions vary with direction, we tested a series of thermal plasma spectra for 
temperatures ranging from 10 5 to 10 6 5 K and searched for a combination of SWCX spectra and thermal 
emission matching the observed intensities and band ratios, while simultaneously being compatible with 
O VI emission measurements. In the frame of collisional equilibrium models and for solar abundances, the 
range we derive for hot gas temperature and emission measure cannot reproduce the Wisconsin C/B band 
ratio. This implies that accounting for SWCX contamination does not remove these known disagreements 
between data and classical hot gas models. We emphasize the need for additional atomic data, describing 
consistently EUV and X-ray photon spectra of the charge-exchange emission of heavier solar wind ions. 



Subject headings: interplanetary medium - ISM: general - ISM: bubbles - supernovae remnants - X-rays: general - 
X-rays: diffuse background - X-rays: ISM 



Introduction 



The diffuse soft X-ray background (SXRB), first 
observed in the 70's (Bowve r et al . 1968; Williamson et al 



1974; Sanders et al. 1977) has since been shown to be 
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the sum of local and distant sources. Above 2 keV it 
is dominated by the extra-galactic background, itself a 
combination of unresolved point sources and warm- 
hot interstellar medi um (WHIM) diffuse emission 
(IHasinger et al.lll993l) . At lo wer energies it is dom- 
inated by the galacti c halo ( Burrows & Mendenhalll 



19911: ISnowden et alJ 119941) . and finally below 0.3 



keV it is mainly due to the unabsorbed emission 
from hot gas filling th e so c a lled Local Bubble 
(LB) dMcCammon et al.l Il983b iBloch et al.l Il986t 



Snowden et al. Il990allbh . a cavity devoid of dense 
gas exten ded at high latitudes and connected to 



the halo dFrisch & Yorkl 1 19831: IWelsh et alJ Il998 



Lallement et al.l 12003b . The main tools used to dis- 
entangle local and distant emission are the 'shadow- 
ing' experiments, i.e. spatial variations of intensity 
and spectral characteristics around and towards dense, 



soft X - ray absorbing clouds (e.g. [ Herbstmeier et al . 
19951) . ISnowdenetS] dl998i l200d) used more than 



370 ROSAT shadows to produce almost full-sky map- 
ping of the 'unabsorbed' component of the emission, 
i.e. the LB contribution. 

This was the generally accepted scenar io until the 



disco very of X-ray emission in comets dLisse et al 



19961) and the identification of the emission mecha 
nism as charge-exchange (CX) reactions between the 
highly ch arged heavy so lar wind ions a nd the cometary 
neutrals dCravendfl997l) . ICoxl dl998l) suggested that 
the CX reactions should also occur between heavy 
SW ions and interstellar neutrals (H and He) in inter- 
planetary space and that the resulting X-ray emission 
(SWCX) should have an impact on the SXRB inter- 
pretation. ICravens estimated that the quiescent 
level of SWCX emission could be of the same order as 
the SXRB component attributed to the LB. 

This interplanetary, heliospheric emission is time- 
dependent because of the intrinsic variable nature of 
the solar wind. Short time-scale variations tend to be 
washed out by integration along the line-of-sight an d 
the size of the emitting region dCravens et alj|200ll) . 
but longer term variations, including those related to 
the solar cycle, can cause more persistent changes in 
the heliospheric emission level. In addition, there is 
a contribution from the Earth's magnetosphere, due to 
charge-transfer with exos pheric neutrals. Such e mis- 
sion, studied in detail by Ro bertson et al. I d2006l) . re- 
acts instantaneously to solar wind variations and mag- 
netosphere shape variations, leading to a high variabil- 
ity and the occurence of high intensity peaks follow- 
ing solar events. The spectral characteristics of some 
spectacular enhancements have been re cently recorded 



by the XMM and Suzaku sa tellites dSnowden et al 
2004 iHenlev & Sheltonll2008l) . For such events both 



magnetospheric and heliospheric contributions may be 
present. 

Very likely most of the sharp increases of terres- 
trial origin have been removed from the ROSAT map 
along with the cleani ng procedure of the L ong Term 



Enhancements (LTEs iSnowden et alJ 119941) . as well 



downwind side of the interstellar flow where the gravi- 
tational cone of focused helium is the most reactive re- 
gion. Indeed, most points on the sky in the ROSAT map 
were observed several times over the course of at least 
two days, allowing identification and removal of peri- 
ods of enhanced emission. A debate is still maintained, 
though, about the actual level of the quasi-stationnary 
heliospheric contribution to the ROSAT maps of unab- 
sorbed emission. This contribution is extremely diffi- 
cult to detect from time variations. On the other hand, 
SCWX and hot gas thermal emission have different 
spectral properties, i.e. the observed spectral informa- 
tion should help to disentangle the two processes. This 
is the subject of the present study. For a recent review 



of all types of SWCX phenomena see iBhardwai et al 
d2007l) . 



The first estimates of the stationary helios pheric 
contribution dCravens! uOOOt iLallementl 120041) were 
based on simplifying assumptions about the spectral 
characteristics of the SWCX emission. Since then the 
existence of the SWCX phenomenon has motivated 
theoretical work on exact ph oton yield values for the 
charge transfer collisio ns dKharchenko & Dalgarnol 
l200Ct|Pepino et alj|2004 . as well as a number of labo- 
ratory experiments devoted to the CX emission mech - 
anism. For a recent review see lWargelin et alJ d2008l) . 

The spatial distribution of magnetosphe ric and he- 



liospheric SWCX emission was modeled bv lRobertson & Cravens 
d2003l) . revealing significant variations in brightness as 
a function of eart h location, line-of-s i ght dir ection, and 



activity phase. Koutroumpa et al. ( 20061) computed 



as some heliospheric increases, especially towards the 



similar maps for a few speci fic energy bands after in- 
cluding IPepino et al.l d2004l) 's det ailed CX emissio n 
spectra for C, N, O, and Ne ions. lLallemenJ d2004l) . 
taking into account the specific viewing geometry of 
ROSAT showed that the heliospheric background in 
the i keV band was nearly isotropic and could have 
been unnoticed in the All-Sky Survey maps, while ac- 
counting for a large portion of the signal, and possibley 
the major part at low galactic latitudes. 

Usi ng both the stationnary and time-dependent 
models iKoutroumpa et al.l d2007l) modelled four high- 
latitude shadowing observations and showed that in the 
3/4 keV band, where the oxygen lines (O VII triplet at 
0.57 keV and O VIII line at 0.65 keV) are dominant, 
the SWCX emission from the heliosphere can account 
for all the unabsorbed, local component of the SXRB, 
with no need of a LB emission. In parallel, the so- 
lar wind contribution to the background and its vari- 
ability have been shown to be responsible for some 
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discrepant measurements dSmith et alj 120071) and for 
supposedly low-energy counterp arts of distant objects 
(lire gman & Lloyd-Dav ies 2006). 

A 10 6 K plasma, however, has very little emission 
in the 3/4 k eV band and main l y emit s in the | keV 
band. The iKoutroumpa et alj ( 120071) results, while 
not requiring any LB emission, therefore do not pre- 
clude the existence of 10 6 K gas. Exact calculations of 
the SWCX spectra and intensities below 0.3 keV are 
mandatory if one wants to disentangle LB hot gas dif- 
fuse emission from the SWCX background. In this pa- 
per we examine the SWCX contribution to the \ keV 
spectral region, compare this contribution to observa- 
tions and also to contributions from hot gas at different 
temperatures. 

Independently of the SWCX contribution, a number 
of results have somewhat contradicted the interpreta- 
tion of the unabsorbed soft X-ray background as the 
LB 10 6 K gas emission. 

i) Data from the NASA EUVE satellite and from 
the dedicated CHIPS mission did not detect the 
EUV emission expected from surroundin g 10 6 K gas 
dJelinskv et al.lll995tlHurwitz et al.ll2005l) . It has been 
suggested that a very low metal abundance may be 
responsible for this non-detection, but the required de- 
pletion level corresponds to the physical state of very 
dense clouds, which is unlikely for 10 6 K, tenuous gas. 

ii) The pressure of this hot gas derived from the X- 
ray background is far above the pressure within the lo- 
cal i nterstellar cloud and other clou ds embedded in the 
LB ( lLallemendll998Hjenkins|[2002l) . 

iii) Low latitude absorption measurements of highly 
charged ions such as Si IV C IV and O VI formed in 
conductive interfaces between the hot (10 6 K) gas and 
embedded cold:warm clouds do not seem to corre- 



than what is predicted by thermal emission mod- 
els. This has been particula rly well demonstrated by 
Bellm & Vaillancourd d2005l) who have made a global 



spond to expectations from the m odels (ISlavin & Frisch 
20021: llndebetouw & Shulll2004 . Column densities of 
Si IV and C IV are too small and line-widths too nar- 



row (IWelsh & Lallementll2005l) . and O VI is detected 
only at the periphery of the Local Cavity, while one 
would also expe ct interfaces between the h ot gas and 
the local clouds (Welsh & Lallemendl2008l) . 

iv) Fundamental discrepancies arise also when 
comparing the Wisconsin sounding rocket survey 
data in the B and C bands, and the ROSAT All- 
Sky Survey data in the Rl and R2 bands. The four 
bands are pictured in figure [[jup per panel. In the 
low energy (0.1-0.2 keV) B band dBloch et al.lll986t 
Snowden et al.lfl994 ), the intensity seems to be higher 



study over the 0.1-0.3keV interval. According to this 
work, a best fit to all energy bands is provided by 
very low metallicity gas at 10 5 85 K, but inspection of 
their results (see their Figure 4) reveals significant dis- 
crepancies between measured and observed ratios for 
this best fit solution. Especially, the B/R12 band ratio 
favours a low temperature -10 5 8 K (the B band inten- 
sity is high and favours a shift of the spectrum towards 
low energies), while the R2/R1 ratio favours temper- 
atures above 10 6 K (R2 is relatively high, favouring a 
shift towards high energies). 

Whether or not the existence of the SWCX back- 
ground can help to explain part or all these contra- 
dictions is a question that has now to be addressed. 
This work is a first step in this direction. In section 
|2] we describe the the SWCX emission and spectral 
model we have developed and how we make use in our 
analysis of the Raymond & Smith (R-S) hot plasma 
model. In section [3] we compute the expected SWCX 
emission and the contribution in ROSAT Rl and R2 
ba nds for each of the 37 8 shadow regions observed 
by ISnowden etall d2000h . We compare the SWCX 
R1+R2 inten sity with the unabsorb ed component de- 
rived by the ISnowden et al. d2000l) shadow analysis 
and discuss the distribution of the discrepancies be- 
tween data and the SWCX model. In section @] we 
compute the SWCX model R2/R1 and B/C band ra- 
tios, as well as the corresponding ratios for hot gas (R- 
S model) in collisional equilibrium within a large tem- 
perature range. We compare the modeled ratios with 
the observed band ratios during the two (Wisconsin 
and ROSAT) surveys. In section|5]we search for a com- 
bination of SWCX and hot gas emission compatible 
with the observed intensities and band ratios and we 
compare those solutions with observational constraints 
from O VI and EUV background measurements. In 
section [6] we discuss the results and draw some con- 
clusions. 

2. Model description 
2.1. SWCX Model 

The basic model calculating the SWCX emission 
in the inner heliosphere was thoroughly presented in 
Koutroumpa et al. (2006); parameters appropriate for 
the ROSAT All-Sky Survey are discussed in [3] We 
calculate self-consistently the neutral H and He den- 
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sity distributions in the inner heliosphere (up to ~100 
AU), in response to solar gravity, radiation pressure 
and anisotropic ionization processes for the two neu- 
tral species. Ionization of H atoms is mainly due to 
their charge-exchange collisions with solar wind (SW) 
protons and He atoms are mostly ionized by solar EUV 
photons and electron impact. We also consider the im- 
pact of CX on the SW ion distributions. This interac- 
tion is described in the following reaction: 

X Q+ + [H,He]^X* (Q - l)+ + [H + ,He + ] (1) 

The collision rate per volume unit R^e+ (cm -3 s ) 
of X G+ ions with the neutral heliospheric atoms is 
given by the equation: 

%e+W = N K Q+(r)v,-{(T (HX Q +) n H (r)+ <T {HeX Q t) n He (r)}2) 

where cr (Hxe+) and <x (He x e +) are the hydrogen and he- 
lium CX cross-sections, n#(r) and n// e (r) are the hy- 
drogen and helium density distributions respectively, 

v T - Vsw — v„ ~ Vsw the relative velocity between 
SW ions and IS neutrals in the inner heliosphere, and 
NxQ+(r) is the self-consistent solution to the differential 
equation: 

dN x a* 



dx 



■ N XQ+ (°"(tf,Xe+) »H(X) + <T(He,xQ+) »»W) 



(3) 



+ Mf<e+n+ (°"(H,A-<e+i>+) n H (x) + a- ( // c> x<0+i)+) nne(x)) 

expressing the evolution of the density distribution of 
ion X e+ along S W streamlines due to production (from 
CX reactions of ion X (e+1)+ ) and loss terms. 

Cross-section uncertainties are mainly due to in- 
strumental systematic errors and most important to the 
collision energy dependance of cross-sections. De- 
tailed uncertainties for individual ions are not given in 
literature, bu t average uncertainti es of ~30% at most 
are reported (IWargelin et al 1 l2008h . 

Then, we establish emissivity grids in units of (pho- 



tons cm 



-3 o-iv 



£;0) = R{XQ+,H)(r) Y(E„H) + R (XO*,He)(r) Y(E h He) (4) 

where Y^M) i s the photon emission yield (in number 
of photons) computed for a spectral line of photon en- 
ergy E, following CX with the corresponding neutral 
species M (H or He individually). For any line of sight 
(LOS) and observation date, the directional intensity 
of this spectral line is given by: 



I Ei (LU) 



4tt X 



■100AU 



Si(s) ds 



(5) 



which defines the average intensity, in Line Units (LU 
= photons cirT 2 s _1 sr _1 ), of the spectral line for the 



particular date and LOS, as well as the solar cycle 
phase (minimum or maximum) corresponding at this 
date. The intensity is somewhat underestimated be- 
cause of the SW ion propagation in the heliosheath up 
to the heliopause, and in the heliotail up to ~3 000 AU, 
where all ions are used up. The outer heliospheric re- 
gion is neglected in our model, but estimates yield a 
maximum additional ~20% contribution in the down- 
wind direction, with possible effects on the SWCX 
spectral hardness (see @. 

Our original atomic database ( Kharchenkol 2005 ) 
included C 5 ' 6+ , N 5A7+ , 6 - 7 - 8+ , Ne 8 - 9+ and Mg 10 - 11+ 
ions. Exact calculations of the cascading photon spec- 
tra were performed individually for these ions when 
they charge exchange with hydrogen and helium re- 
spectively. Detailed CX collision cross sections taking 
into account both the neutral target species and the 
solar wind velocity regime were include in the cal- 
culations (P. Stancil private communication). These 
calculations have already been used to reproduce 
observed SWCX sp ectra from comets with CHIPS 
dSasseen et al.ll2006l) . 

The database was recently updated to include 
p e 7...i3 +) si 5 - 10+ , S 6 ~ 11+ , Mg 4 - 9+ ions that emit in- 
tense lines in the 0.1-0.3 keV range. Individual emis- 
sion spectra induced in the charge exchange collisions 
of these ions have very complicated structures be- 
cause of a large number of intermediate multiplets 
related to different excited states of many-electron 
ions. The photon yields Y(E,M) for heavier ions 
were calculated using the simpli fied hydrogenic model 
( Kharchenko & Dalgarnoll2000l) . which assumes a hy- 
drogenic nature of electronic excited states. In this 
model, the effective charge of hydrogenic ion is com- 
puted from an accurate value of the ion ionization 
potential and branching ratios of radiative cascading 
transitions are chosen to be the same as in all H-like 
ions. Moreover, photon yields were calculated using 
a single neutral species, which means that no distinc- 
tion between H and He was made. The hydrogenic 
approximation of the CX emission spectra is a quan- 
tum mechanical model in which an actual ion spec- 
tra may be replaced with the hydrogenic spectra. In 
this model the total energy of emitted photons is de- 
fined by an initial state-population and should be an 
accurate quantity matching real spectra. Positions of 
emission lines do not correspond exactly to real emis- 
sion spectra, but this defect is not very important at 
the low resolution of the observed spectra. Total cross 
sections of CX collisions for the hydrogenic approx- 
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Fig. 1. — Top panel: Wisconsin B (plain) and C (dashed), and ROSAT Rl (dash-dotted) and R2 (dotted) band effective 
areas. Bottom panel: Example of calculated SWCX spectra in Line Units (photons cirT 2 s _1 sr _1 ). The emitting ions 
are marked above the most prominent lines. 



imation s have been calculated using t he over-barrier 
model dKharchenko & Dalgarnoll200lb . An example 
of calculated spectra is presented in figure flower 
panel, with the emitting ion identifying the most in- 
tense lines. 

2.2. Hot Gas thermal emission 

We use a Raymond-Sm ith (RS) hot plasma model 
dRaymond & Smith! 1 19771) assuming typical metal 
abundances [He, C, N, O, Ne, Mg, Si, S, Ar, Ca, 
Fe, Ni] = [10.93, 8.52, 7.96, 8.82, 7.96, 7.52, 7.60, 
7.20, 6.90, 6.30, 7.60, 6.30] dAllerJl973h. We use this 



model rather than the APEC model dSmifh et alj|2001 



that superseded it because we are most interested in 
the j keV range. APEC includes only transitions for 
which accurate atomic rates are available, while the 
code of Raymond and Smith estimates the emission in 
the large number of weak lines that are known to be 
present (e.g., from moderately ionized species of Mg, 
Si, S and Fe) but which lack accurate excitation rates 



and wavelengths. Given the low spectral resolution of 
the observations considered here and our interest in the 
total emitted power, the RS model serves very well. 

This model gives us X-ray emissivities fi(T) and 
f2(T) convolved by and summed in the ROSATRl and 
R2 bands, respectively, as a function of temperature 
such that the total hot gas X-ray flux in these bands is 
defined as: 

hl,LB = h,LB + h,LB = EM(i, T) ■ (fl(T) + f 2 (T)) (6) 

where EM(i, T) is the emission measure for tempera- 
ture T and look direction i. Units of functions fi and 
fo are RU EM -1 , where RU = 10~ 6 cts s _1 arcmin~ 2 is 
the usual ROSAT detector unit and EM is the typical 
emission measure unit cirT 6 pc. 

Equivalently, the hot gas emission in bands B and 
C is defined as: 

hc,LB = Ib,lb + Ic,lb = EM(i,T) ■ (f B (T) + f c (T)) (7) 
where ffl(T) and fc(T) are the equivalent emissivity 
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ABSOLUTE GALACTIC LATITUDE (°) 



Fig. 2. — Measured unabsorbed (112,0/™, from lSnowden et al. (2000), gray circles) and calculated SWCX (black dots) 
fluxes in the ROSAT j keV (R1+R2) band as a function of absolute galactic latitude for 378 shadowing observations. 
Flux is given in ROSAT Units (10~ 6 cts s arcmkr 2 ). 



functions in the B and C bands derived by the RS 
plasma code, in units of cts s EM -1 . 

3. R1+R2 intensities 

We have calculated SWCX spectra in ROSAT ob- 
servation geometry for th e shadow field l i nes of sight 
(LOS) listed in table 1 of ISnowden et al.l j2000h . that 
were observed during the ROSAT all-sky su rvey. The 
shadows analysed by ISnowden et al.l (|2000) were lo- 
cated at high galactic latitudes and in general above 
15° from the galactic plane. The ROSAT observation 
geometry is defined with the view direction perpen- 
dicular to the Sun-satellite direction. Thus, it takes a 
six-month period to build a full-sky map of the soft 
X-ray intensity. The ROSAT all-sky survey was per- 
formed between July 1990 and February 1991, which 
corresponds to maximum solar activity conditions that 
were taken into account in the SWCX simulations. 

Maximum solar activity conditions imply the fol- 
lowing input parameters in the SWCX model. We 
consider a radiation pressure to gravity ratio jj = 1.5 
for neutral hydrogen and slightly anisotropic ioniza- 
tion rates varying between 8.4xl0~ 7 s -1 at the solar 



verse ionization rate) at 1 A U is 0.62xl0 7 s at sola r 
maximum, in agreement with iMcMullin et al.l (120041) . 
In solar maximum, the solar wind is considered to be 
a complex mix of slow and fast wind states that is in 
general approximated with an average slow wind flux. 
Slow solar wind flows at ~ 400 km/s and has a proton 
density of ~ 6.5 crrT 3 at the Earth position. The oxy- 
gen content with respect to protons is [O/H] = 1/1780. 
The most important heavy ion charge state abundances 
with respect to oxygen [X q+ / O] are: C 5M : [0.21, 
0.318], 6 ' 7 ' 8+ : [0.73, 0.2, 0.07], Si 8A10+ : [0.057, 
0.049, 0.021] and Fe 8A ' a " + : [0.034, 0.041, 0.0 31, 
0.023] (adopted from ISchwadron & Cravensll2000l) . 



equator and 6.7x10 s at the poles dOuemerais et al 



2006). For neutral helium, the average lifetime (in- 



We have convolved the individual spectra with the 
ROSAT Rl and R2 band responses in order to calculate 
the total SWCX flux in these bands, as well as the total 
R12 (R1+R2) flux. We plot the resulting R12 SWCX 
flux and the unabsor bed 1 12,0/™ component from the 
Sno wden et al. I (120001) analysis as a function of abso- 
lute galactic latitude in figure[2] Ii2,o»j corresponds to 
the unabsorbed portion of the SXRB that was origi- 
nally attributed to the LB ~10 6 K hot gas. X-ray in- 
tensities are presented in ROSAT Units (RU = 10~ 6 cts 
s _1 arcmin~ 2 ). 

The SWCX R12 flux (black dots) varies between 
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Fig. 3. — Unabsorbed ln,obs R12 flux data over simulated lu,cx SWCX flux ratio map in galactic coordinates with an 
equal-area Aitoff projection. For information we note the UW and DW directions in the sky. The LC chimneys appear 
clearly where the In.obs/hi.cx ratio is the highest. 



212 and 460 RU with an average value of 332 RU and 
is fairly uniformly distributed across all latitudes. The 
lower and upper limits calculated in the SWCX sim- 
ulations for average maximum conditions are repre- 
sented with the plain black horizontal lines. 

On the other hand, the una bsorbed T ?„/,., compo- 
nent (gray circles) derived in the lSnowden et al. I d2000h 
analysis has a clear correlation with the absolute galac- 
tic latitude. Higher ln,obs values are measured towards 
higher latitudes, where the local cavity is enlarged and 
communicates with the galactic halo through the chim- 
neys. 

In the figure it is clear that the SWCX intensity is 
of the same order as the l\2,obs intensity measured in 
low galactic latitudes (up to around 20-25°). We can 
conclude, then, that the SWCX 4 keV flux could ac- 
count for most of the observed ROSAT emission in the 
galactic plane. This conclusion implicitly assumes that 
the highly peaked exospheric SWCX contribution has 
been cleaned from the ROSAT data, but not the helio- 
spheric contribution. 



Fig[3]is a map in galactic coordinates of the ratio be- 
tween the unabsorbed ROSAT emission and the com- 
puted SWCX contribution. The map clearly reveals 
the emission from the so-called chimneys that connect 
the Local Cavity to the northern and the southern halo. 
Our intensity results do not preclude that outside these 
chimneys the totality of the signal is SWCX emission. 

4. Band ratios 

For each SWCX spectrum calculated in the look di- 
rections presented in figure [2] we have calculated the 
R2/R1 (ROSAT) and C/B0 (Wisconsin) ratios. We find 
an average R2/R1 (hereafter Rex) ratio of 1.39 and an 
average C/B ratio of 6.67. Although both ratios show 
very little variation across the sky, there is a hardness 
trend of the SWCX spectra with harder spectra towards 
the downwind direction (UW to DW variations: Rex 



1 The original papers on the Wisconsin survey referred to the B/C ra- 
tio, but given the extensive use of the R2/R1 ratio in our analysis 
and the rough correspondence between Rl and B, and R2 and C, we 
refer to the C/B ratio. 
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Fig. 4. — Wisconsin (upper panel) and ROSAT (lower 
panel) band ratios for data (UW-bold, DW-plain line), 
SWCX model (dashed line) and the R-S hot gas model 
(dotted line) as a function of temperature. 



= [1.36 - 1.41], C/B = [6.25 - 7.14]). However, we 
need to alert the reader that these are somewhat uncer- 
tain SWCX spectra and therefore somewhat uncertain 
band ratios. Indeed, as we mentioned in section|2l ex- 
act calculations for Fe, Si, Mg, Al are not yet avail- 
able, and no distinction was done between the neutral 
targets (H or He), while laboratory experiments show 
that the energy levels populated after the electron cap- 
ture and the subsequent radiative cascades may differ 
significantly for different targets. Since most of the 
SWCX downwind emission is due to the interaction 
with neutral helium, while on the upwind side hydro- 
gen is the main contributor, more precise calculations 
could have an effect on the hardness. Also, although 
preliminary calculations show an almost negligible ef- 
fect, the heliospheric model cutoff (especially in the 
DW directions) may be responsible for the "loss" of 
relatively more emission from lower charge states (at 



relatively lower energies) than emission from higher 
charge states (at relatively higher energies). Thus, the 
calculated SWCX spectra may actually be softer than 
what is predicted here. It is evident that a more de- 
tailed calculation taking into account all metals and 
the neutral target nature, as well as detailed cascad- 
ing collisions (secondary ion production) in the outer 
heliosphere is needed in the future. On the other hand, 
the interval we find for the ratio can be used a reliable 
value for the average SWCX ratio. 

These SWCX ratios have to be compared with the 
corresponding ratios for thermal emission. The latter 
were obtained as a function of temperature by con- 
volving the Raymond-Smith spectra with the ROSAT 
band responses Rl and R2 and Wisconsin B and C re- 
sponses. The results are shown in figure [4] Above 
logT= 6.1 the thermal R2/R1 ratio reaches its maximal 
value of 1.2. It remains however slightly lower than 
the SWCX ratio of 1.36-1.41. At those temperatures 
the thermal C/B ratio increases to its maximal value of 
=i 4, a value almost half the SWCX ratio of 6.24-7.14, 
i.e. a significant difference. Those curves allow esti- 
mates of the ratios for combinations of thermal plus 
SWCX background emissions. 

5. Combination of the heliospheric SWCX and 
LB hot plasma emission 

For a comparison with the data we consider two re- 
gions: one centered on the direction of the incoming 
IS flow at ecliptic coordinate s (A, P) ~ (252.3 8.5)° 



for the IS H flow, according to lLallement et alj d2005l) 
(upwind -UW- direction) and one looking at the out- 
going flow direction (downwind -DW- direction). In 
galactic coordinates the UW direction corresponds to 
(1, b) = (5.4, 18.9)°, close to the galactic center di- 
rection (anti-galactic direction for DW respectively). 
These two regions are also very close to the minimum 
an d maximum va l ues of the hardness ratio derived 
by lSnowden et al. I dl990bh . which define the so-called 
color gradient axis of the soft X-ray background. For 
these two regions we can derive average values of ob- 
served unabsorbed j k eV emission us i ng th e Local 



Bubble contours in the ISnowden et al.l d 19981) analy- 
sis. For the upwind (UW) direction the ROSAT un- 
absorbed In.obs emission we estimate ~325 RU, while 
for the downwind (DW) direction the observed unab- 
sorbed level is found to be 450 RU. 

The equivalent B+C intensities in the Wiscon- 
sin survey are estimated on average ~90 cts s _1 and 
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~125 c ts s 1 for the UW and DW directions respec- 



tively ( Snowden et al. 1990bl) . However, those in- 



tensities include both the foreground (assumed LB) 
and more distant components (galactic halo and ex- 
tragalactic), since the Wisconsin survey did not have 
enough spatial resolution to study the shadowing 
fields. Moreover, the Wisconsin sounding rocket mea- 
surements looked in the roughly anti- Sunward direc- 
tion, which should also affect the comparison with the 
ROSAT All-Sky Survey in terms of the SWCX compo- 
nent spatial distribution. For instance, for the DW look 
directions, the Wisconsin sounding rockets were ob- 
serving directly through the He cone and should have 
had a higher "contamination" of SWCX emission than 
the ROSAT detectors that must have been located in 
crosswind positions on the Earth's orbit in order to 
observe in the DW directions. 

The measured R2/R1 and C/B ratios for these re- 
gions are shown superimposed to the models in figure 
[4] and listed in Table 1. For the UW area the measured 
R2/R1 ratio is close to the SWCX value, although 
slightly smaller, while for the DW area it is signifi- 
cantly smaller. For both areas the measured C/B ratio 
is lower than the SWCX ratio. Fig |4] shows that for 
both areas a combination of SWCX and thermal emis- 
sions may in principle account for the ROSAT mea- 
surements, the thermal emission lowering the R2/R1 
ratio to achieve the observed value. Similarly, inde- 
pendently of the R2/R1 ratio, fig [4] also shows that a 
combination of both backgrounds may account for the 
Wisconsin data, the thermal emission decreasing the 
C/B ratio to achieve the observed value. 

It remains to find a combination satisfying both ra- 
tios simultaneously. Our attempt to find a solution is 
the following one. For each assumed temperature of 
the hot gas, we use the Rl and R2 data (and thus the 
observed ratio) and the SWCX spectral shape to de- 
rive the respective contributions of SWCX and thermal 
emission, i.e., we derive which quantity of SWCX in- 
duced R12 intensity and which emission measure EM 
for the hot gas lead to the measured intensities and the 
measured R2/R1 ratio. We a posteriori calculate the B 
and C intensities, C/B ratio and the O VI emission of 
the hot gas and compare with the data. 

The SWCX model predicts a total R12 intensity: 



'12.CX 



= / 



i,cx + h,cx 



= h.cx ■ (1 + Rex) (8) 
1 +Rcx 



h,cx 



Rex 



(9) 



model. 

The total unabsorbed flux l\2,obs& measured in the 
R12 band towards look direction i is the sum of LB hot 
gas Ii2,lb and SWCX I 12 , C x fluxes: Ii 2 , & s (i) = hi.cx + 
h2,LB, so that SWCX intensity can be written: 

hi,cx = hiMi)-hi,LB = In,obs(i)-EM-(fi (T)+f 2 (T)) 

(10) 

The observed R2/R1 ratio (hereafter Robs) towards 
look direction i is defined by the equation: 



R o bs(0 



h,cx + h.LB 
h,cx + h,LB 



(11) 



Resolving equation QT] by using equations l6l to [TUl 
we find the hot gas emission measure EM(i, T) as a 
function of temperature, total R12 measured intensity 
Ii2,oi«(i) an d measured R„/ M (i) ratio towards look direc- 
tion i. 



EM(i, T) 



Rex ~ Robs(i) 



{R obs (i) + 1) • (Rcx-fdT) - fz(T)) 



'12,0/7, 



.(/) 



where R cx is the R2/R1 ratio predicted by the SWCX 



(12) 

We calculate the emission measure for the two UW and 
DW directions defined above using the following nu- 
merical values: (i) the Rex ratio is constant and equal 
to 1.39, (ii) observed values of the unabsorbed portion 
of the i keV emission in the R12 band are Ii2,oi. s (UW, 
DW) = (3 25, 450) RU as derive d from the LB con- 
tours in the lSnowden et alJd 19981) analysis for the UW 
and DW (respectively galactic and anti-galactic) direc- 
tions, and (iii) the corresponding observed R2/R1 ratio 
is R obs = 1.25 and 1.04 for the UW and DW directions 
respectively. 

In figure [5] we show the resulting EM and the por- 
tion of the total emission due to the SWCX mechanism 
(called residual emission) as derived from equations 
[T2l and [TO] respectively as a function of logT and for 
the two look lines. We show superimposed the R12 
measured intensities in those directions. For the cal- 
culated EM and corresponding temperatures we have 
also added to figure|5]the intensity of the O VI doublet 
at A = 1034 A (1032 A and 1038 A). In order to cal- 
culate this O VI doublet emissio n we have used equa- 
tion (5) of lShull & Slavinl d 1994 and assumed that in- 
terstellar O abundance is 8.5x 10~ 4 . We also assume 
that the O VI ion proportion depends on temperature 
according to the C hianti database fo rmulae for colli- 
sional equilibrium dLandi et alJuOQq) . 

In order to delimit the possible temperature solu- 
tions for the LB hot gas, we place the following con- 
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Fig. 5. — Solutions for the EM (dashed line) and residual SWCX emission (dot-dashed line) as a function of tem- 
perature of a hot plasma for observed R12 intensities (Li2,obs> plain horizontal lines) for an upwind (left panel) and 
downwind (right panel) look directions. Upper O VI doublet intensity limit (dotted line) is calculated for the EM 
values. Discontinuities in the O VI intensity curve are due to O VI density interpolation i n the Chianti d atabase. The 
cross-hatched regions show the hot gas temperature ranges excluded by O VI intensity JSheltonll2003l) and SWCX 
emission constraints. 



straints: i) We assume that the SWCX model is accu- 
rate enough to ensure that the heliospheric emission in 
the R12 band cannot be lower than ~ 212 RU (lower 
limit in figure [2]). This gives us (from the right panel 
of FigO an upper limit of logT = 6 in temperature, ii) 
We use the observed upper li mit of O VI dou blet in- 
tensity, reported at ~ 800 LU dSheltonll2003h . which 
gives us two temperature limits at logT = 5.13 and 
5.64 (extreme limits in the DW direction). The inter- 
val between those two temperatures is forbidden be- 
cause the corresponding O VI column densities (and 
thus the O VI intensity) are too high to match observa- 
tions. Temperatures below 10 513 K would predict ex- 
tremely strong C VI and N VI absorption t oward nearby 



stars, which have not been obse rved (e.g. lLehner et al. 
2003l:IWelsh & Lallementl2"005l) so we do not consider 
it a realistic solution. The limits of valid temperature 
intervals are marked by the vertical bold lines and the 
cross-hatched regions in figure [5] show the excluded 
temperature ranges. The two most plausible hot gas 
temperature limits (logT = 5.64, 6.0) along with the 
corresponding UW and DW emission measures and 
residual SWCX emission are summarized in table [1] 

We also calculate for the two boundary solutions 
the corresponding SWCX intensities and the thermal 
emission intensities in the B and C bands by convolv- 
ing our simulated SWCX spectra and the hot gas spec- 
tra with the band responses. For the two temperatures 
the total hot gas and SWCX intensity in (B + C) band 



is found to be about 83 cts s and about 1 15 cts s for 
the UW and DW directions respectively. This similar- 
ity arises from the similarity between the wavelength 
intervals covered by the B and C bands and the Rl and 
R2 bands (see fig. [TJ. The C/B ratio derived from this 
analysis is 5.0 and ~3.5 for the UW and DW directions 
accordingly, for temperatures above 10 5 54 K. 

The (B+C) total intensity is consistent with the 

lower v alues reported in the Wisconsin survey dSnowden et al 
1990bh . which c orrespond to the lower galactic lati 



tudes. Moreover, Snowden et al. ( 1990bl) did not pro- 
ceed with a shadowing analysis of the Wisconsin data, 
so the reported values include both local and more 
distant absorbed components and are expected to be 
higher than the hot gas and SWCX combination we 
present here. 

However, the C/B ratio computed in the analysis (5 
to ~3.6 from UW to DW, depending on temperature) 
is inconsistent with the observed value, especially in 
the DW direction (observed ~2.2), suggesting that we 
should need more hot gas emitting in the B band. This 
inconsistency cannot be attributed to the absorbed por- 
tion of emission included in the Wisconsin data anal- 
ysis because the absorbed component is a high-T gas 
giving a harder spectrum since absorption is more ef- 
fective in lower energies. 

This inconsistency of the DW C/B ratio is impor- 
tant, since it seems difficult to explain in the context 
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Table 1 : Temperature (logT) limits and E.M. solutions for upwind (UW) and downwind (DW) look directions when 
combining a SWCX and RS hot plasma code. In the lower part of the table we include the observational input 
considered. 
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of our study. As a matter of fact, as can be seen in 
table [1] the SWCX contribution in the C band is large 
whatever the temperature, and reaching a C/B ratio of 
2.2 requires a very small SWCX emission, in our sense 
far from realistic. Again, as we discussed in the intro- 



duction and it was shown in the Bellm & Vaillancourt 



(12005b study, the B intensity is higher than expected 
from the models. This seems to remain true (and even 
worse) when taking into account the SWCX contribu- 
tion. 

6. Discussion and conclusions 

We have modeled the intensity and spectral char- 
acteristics of the heliospheric SWCX emission at the 
time of the ROSAT survey and compar ed with the un- 
absorbed, local emission derived by ISnowden et al 



d2000h in the j keV band. The results show that 
the SWCX emission can account for most of the to- 
tal intensity recorded in the R1+R2 bands for most of 
low latitude lines-of-sight. A map of the heliospheric 
SWCX portion of the total signal clearly reveals the 
high latitude chimneys to the halo as the only regions 
unambiguously dominated by hot gas emission. Such 
a result can be interpreted as meaning that little or no 
hot gas exists within the galactic disk. 



The spectral characteristics however reveal more 
complexity and preclude such a simple scenario. The 
SWCX band ratios disagree with the observations, es- 
pecially towards the galactic anti-center and at low en- 
ergies (C/B). We have thus searched for a combination 
of SWCX and thermal emission from hot gas in equi- 
librium and solar abundances able to reproduce the 
data. Our study shows that a combination of SWCX 
and thermal emission can reproduce the data in the 
galactic center hemisphere at low latitudes. For this 
solution the SWCX emission strongly dominates. The 
temperature of the hot gas is constrained within the 
interval 10 5 64 -10 6 . The upper limit is constrained by 
the lower limit on the SWCX intensity. This upper 
limit can be considered as firmly determined, thanks 
to recent observational studies above 0.3 keV that have 



confirmed the validity of our model (IKoutroumpa et al 



20071) . The temperatures lower than 10 5 64 are ex 



eluded by O VI emission observations ( Sh eltonll2003l) 
and interstellar ion absorption lines toward nearby 
stars. 

On the other hand, it is difficult to fit with such a 
combination the Wisconsin data. In the UW (galactic 
center direction) a combination of hot gas and SWCX 
emission gives (B+C) intensities as well as C/B ratios 
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roughly compatible with the observed values for sev- 
eral different temperature ranges. The main difficulty 
is the impossibility to account for the very low C/B 
ratio measured towards the galactic anti-center direc- 
tion with the present input models used in our study. 
We note that the high B intensity is also clearly a 
problem for any hot gas solution, including the very 
high depletion hypothesis, a s shown by the study of 
Bellm & Vaillancourtl (12005b . The SWCX contribu- 
tion, which hardens the spectra, reinforces this diffi- 
culty. 

However, further investigation is required on the 
model's uncertainties in order to quantify their in- 
fluence on the spectral hardness of SWCX emission. 
Further analysis of the hydrogenic ion approximation 
is needed, since hydrogenic ions tend to emit pho- 
tons at higher energies following CX (since high-n 
to ground transitions are generally allowed) than the 
multi-electron ions considered here (because selection 
rules and more complicated atomic structure lead to 
more cascades before the final transition to ground. 
Therefore, it is likely that the model spectra have sig- 
nificantly less flux at low energies than they should. 
Given the steeply decreasing effective areas at lower 
energies (see Fig. [1), this would help explain some of 
the discrepancies seen in the C/B ratio. Also, the fact 
that no distinction was made between the neutral tar- 
gets (H or He) in the \ keV calculations, would also 
have an effect on the spectral hardness, since there is 
an effect on the electron capture level (roughly pro- 
portional to V13.6eV//,„ where /„ is the neutral tar- 
get ionization potential), the capture level with He be- 
ing somewhat lower than with H. The ROSAT obser- 
vation geometry tends to smooth out those differences, 
but the SWCX spectrum will tend to be softer in the 
DW directions (i.e., looking through or near the He 
cone), again helping to explain some of the C/B ra- 
tio anomaly. Finally, at the low energies of the B 
band, more detailed calculations of the heliospheric 
and magnetospheric signals must be performed, espe- 
cially the low energy secondary SWCX emissions in 
the heliosheath and heliotail, i.e. subsequent recombi- 
nations of partially neutralized solar wind high ions. 
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